Post-transcriptional splicing can occur in a slow-moving zone around the gene.

Splicing is the stepwise molecular process by which introns are removed from pre-mRNA and exons are joined together to form mature mRNA sequences. The ordering and spatial distribution of these steps remain controversial, with opposing models suggesting splicing occurs either during or after transcription. We used single-molecule RNA FISH, expansion microscopy, and live-cell imaging to reveal the spatiotemporal distribution of nascent transcripts in mammalian cells. At super-resolution levels, we found that pre-mRNA formed clouds around the transcription site. These clouds indicate the existence of a transcription-site-proximal zone through which RNA move more slowly than in the nucleoplasm. Full-length pre-mRNA undergo continuous splicing as they move through this zone following transcription, suggesting a model in which splicing can occur post-transcriptionally but still within the proximity of the transcription site, thus seeming co-transcriptional by most assays. These results may unify conflicting reports of co-transcriptional versus post-transcriptional splicing.

On a kneading theory for gene-splicing.

Two well-known facets in protein synthesis in eukaryotic cells are transcription of DNA to pre-RNA in the nucleus and the translation of messenger-RNA (mRNA) to proteins in the cytoplasm. A critical intermediate step is the removal of segments (introns) containing ∼97% of the nucleic-acid sites in pre-RNA and sequential alignment of the retained segments (exons) to form mRNA through a process referred to as splicing. Alternative forms of splicing enrich the proteome while abnormal splicing can enhance the likelihood of a cell developing cancer or other diseases. Mechanisms for splicing and origins of splicing errors are only partially deciphered. Our goal is to determine if rules on splicing can be inferred from data analytics on nucleic-acid sequences. Toward that end, we represent a nucleic-acid site as a point in a plane defined in terms of the anterior and posterior sub-sequences of the site. The "point-set" representation expands analytical approaches, including the use of statistical tools, to characterize genome sequences. It is found that point-sets for exons and introns are visually different, and that the differences can be quantified using a family of generalized moments. We design a machine-learning algorithm that can recognize individual exons or introns with 91% accuracy. Point-set distributions and generalized moments are found to differ between organisms.

Mutations and intron polymorphisms in voltage-gated sodium channel genes of different geographic populations of Culex pipiens pallens/Culex pipiens quinquefasciatus in China.

BACKGROUND: Culex pipiens pallens and Culex pipiens quinquefasciatus are the dominant species of Culex mosquitoes in China and important disease vectors. Long-term use of insecticides can cause mutations in the voltage-gated sodium channel (vgsc) gene of mosquitoes, but little is known about the current status and evolutionary origins of vgsc gene in different geographic populations. Therefore, this study aimed to determine the current status of vgsc genes in Cx. p. pallens and Cx. p. quinquefasciatus in China and to investigate the evolutionary inheritance of neighboring downstream introns of the vgsc gene to determine the impact of insecticides on long-term evolution. METHODS: Sampling was conducted from July to September 2021 in representative habitats of 22 provincial-level administrative divisions in China. Genomic DNA was extracted from 1308 mosquitoes, the IIS6 fragment of the vgsc gene on the nerve cell membrane was amplified using polymerase chain reaction, and the sequence was used to evaluate allele frequency and knockdown resistance (kdr) frequency. MEGA 11 was used to construct neighbor-joining (NJ) tree. PopART was used to build a TCS network. RESULTS: There were 6 alleles and 6 genotypes at the L1014 locus, which included the wild-type alleles TTA/L and CTA/L and the mutant alleles TTT/F, TTC/F, TCT/S and TCA/S. The geographic populations with a kdr frequency less than 20.00% were mainly concentrated in the regions north of 38° N, and the geographic populations with a kdr frequency greater than 80.00% were concentrated in the regions south of 30° N. kdr frequency increased with decreasing latitude. And within the same latitude, the frequency of kdr in large cities is relatively high. Mutations were correlated with the number of introns. The mutant allele TCA/S has only one intron, the mutant allele TTT/F has three introns, and the wild-type allele TTA/L has 17 introns. CONCLUSIONS: Cx. p. pallens and Cx. p. quinquefasciatus have developed resistance to insecticides in most regions of China. The neighboring downstream introns of the vgsc gene gradually decreased to one intron with the mutation of the vgsc gene. Mutations may originate from multiple mutation events rather than from a single origin, and populations lacking mutations may be genetically isolated.

Intron detention tightly regulates the stemness/differentiation switch in the adult neurogenic niche.

The adult mammalian brain retains some capacity to replenish neurons and glia, holding promise for brain regeneration. Thus, understanding the mechanisms controlling adult neural stem cell (NSC) differentiation is crucial. Paradoxically, adult NSCs in the subependymal zone transcribe genes associated with both multipotency maintenance and neural differentiation, but the mechanism that prevents conflicts in fate decisions due to these opposing transcriptional programmes is unknown. Here we describe intron detention as such control mechanism. In NSCs, while multiple mRNAs from stemness genes are spliced and exported to the cytoplasm, transcripts from differentiation genes remain unspliced and detained in the nucleus, and the opposite is true under neural differentiation conditions. We also show that m6A methylation is the mechanism that releases intron detention and triggers nuclear export, enabling rapid and synchronized responses. m6A RNA methylation operates as an on/off switch for transcripts with antagonistic functions, tightly controlling the timing of NSCs commitment to differentiation.

PPR596 Is Required for nad2 Intron Splicing and Complex I Biogenesis in Arabidopsis.

Mitochondria are essential organelles that generate energy via oxidative phosphorylation. Plant mitochondrial genome encodes some of the respiratory complex subunits, and these transcripts require accurate processing, including C-to-U RNA editing and intron splicing. Pentatricopeptide repeats (PPR) proteins are involved in various organellar RNA processing events. PPR596, a P-type PPR protein, was previously identified to function in the C-to-U editing of mitochondrial rps3 transcripts in Arabidopsis. Here, we demonstrate that PPR596 functions in the cis-splicing of nad2 intron 3 in mitochondria. Loss of the PPR596 function affects the editing at rps3eU1344SS, impairs nad2 intron 3 splicing and reduces the mitochondrial complex I's assembly and activity, while inducing alternative oxidase (AOX) gene expression. This defect in nad2 intron splicing provides a plausible explanation for the slow growth of the ppr595 mutants. Although a few P-type PPR proteins are involved in RNA C-to-U editing, our results suggest that the primary function of PPR596 is intron splicing.

Development of Genome-Wide Intron Length Polymorphism (ILP) Markers in Tea Plant (Camellia sinensis) and Related Applications for Genetics Research.

The market value of tea is largely dependent on the tea species and cultivar. Therefore, it is important to develop efficient molecular markers covering the entire tea genome that can be used for the identification of tea varieties, marker-assisted breeding, and mapping important quantitative trait loci for beneficial traits. In this study, genome-wide molecular markers based on intron length polymorphism (ILP) were developed for tea trees. A total of 479, 1393, and 1342 tea ILP markers were identified using the PCR method in silico from the 'Shuchazao' scaffold genome, the chromosome-level genome of 'Longjing 43', and the ancient tea DASZ chromosome-level genome, respectively. A total of 230 tea ILP markers were used to amplify six tea tree species. Among these, 213 pairs of primers successfully characterize products in all six species, with 112 primer pairs exhibiting polymorphism. The polymorphism rate of primer pairs increased with the improvement in reference genome assembly quality level. The cross-species transferability analysis of 35 primer pairs of tea ILP markers showed an average amplification rate of 85.17% through 11 species in 6 families, with high transferability in Camellia reticulata and tobacco. We also used 40 pairs of tea ILP primers to evaluate the genetic diversity and population structure of C. tetracocca with 176 plants from Puan County, Guizhou Province, China. These genome-wide markers will be a valuable resource for genetic diversity analysis, marker-assisted breeding, and variety identification in tea, providing important information for the tea industry.

Genome-Wide Identification of the 14-3-3 Gene Family and Its Involvement in Salt Stress Response through Interaction with NsVP1 in Nitraria sibirica Pall.

14-3-3 proteins are widely distributed in eukaryotic cells and play an important role in plant growth, development, and stress tolerance. This study revealed nine 14-3-3 genes from the genome of Nitraria sibirica Pall., a halophyte with strong salt tolerance. The physicochemical properties, multiple sequence alignment, gene structure and motif analysis, and chromosomal distributions were analyzed, and phylogenetic analysis, cis-regulatory elements analysis, and gene transcription and expression analysis of Ns14-3-3s were conducted. The results revealed that the Ns14-3-3 gene family consists of nine members, which are divided into two groups: ε (four members) and non-ε (five members). These members are acidic hydrophilic proteins. The genes are distributed randomly on chromosomes, and the number of introns varies widely among the two groups. However, all genes have similar conserved domains and three-dimensional protein structures. The main differences are found at the N-terminus and C-terminus. The promoter region of Ns14-3-3s contains multiple cis-acting elements related to light, plant hormones, and abiotic stress responses. Transcriptional profiling and gene expression pattern analysis revealed that Ns14-3-3s were expressed in all tissues, although with varying patterns. Under salt stress conditions, Ns14-3-3 1a, Ns14-3-3 1b, Ns14-3-3 5a, and Ns14-3-3 7a showed significant changes in gene expression. Ns14-3-3 1a expression decreased in all tissues, Ns14-3-3 7a expression decreased by 60% to 71% in roots, and Ns14-3-3 1b expression increased by 209% to 251% in stems. The most significant change was observed in Ns14-3-3 5a, with its expression in stems increasing by 213% to 681%. The yeast two-hybrid experiments demonstrated that Ns14-3-3 5a interacts with NsVP1 (vacuolar H+-pyrophosphatase). This result indicates that Ns14-3-3 5a may respond to salt stress by promoting ionic vacuole compartmentalization in stems and leaves through interactions with NsVP1. In addition, N. sibirica has a high number of stems, allowing it to compartmentalize more ions through its stem and leaf. This may be a contributing factor to its superior salt tolerance compared to other plants.

Highly Reactive Group I Introns Ubiquitous in Pathogenic Fungi.

Systemic fungal infections are a growing public health threat, and yet viable antifungal drug targets are limited as fungi share a similar proteome with humans. However, features of RNA metabolism and the noncoding transcriptomes in fungi are distinctive. For example, fungi harbor highly structured RNA elements that humans lack, such as self-splicing introns within key housekeeping genes in the mitochondria. However, the location and function of these mitochondrial riboregulatory elements has largely eluded characterization. Here we used an RNA-structure-based bioinformatics pipeline to identify the group I introns interrupting key mitochondrial genes in medically relevant fungi, revealing their fixation within a handful of genetic hotspots and their ubiquitous presence across divergent phylogenies of fungi, including all highest priority pathogens such as Candida albicans, Candida auris, Aspergillus fumigatus and Cryptococcus neoformans. We then biochemically characterized two representative introns from C. albicans and C. auris, demonstrating their exceptionally efficient splicing catalysis relative to previously-characterized group I introns. Indeed, the C. albicans mitochondrial intron displays extremely rapid catalytic turnover, even at ambient temperatures and physiological magnesium ion concentrations. Our results unmask a significant new set of players in the RNA metabolism of pathogenic fungi, suggesting a promising new type of antifungal drug target.

The splicing regulators RBM5 and RBM10 are subunits of the U2 snRNP engaged with intron branch sites on chromatin.

Understanding the mechanisms of pre-mRNA splicing is limited by the technical challenges to examining spliceosomes in vivo. Here, we report the isolation of RNP complexes derived from precatalytic A or B-like spliceosomes solubilized from the chromatin pellet of mammalian cell nuclei. We found that these complexes contain U2 snRNP proteins and a portion of the U2 snRNA bound with protected RNA fragments that precisely map to intronic branch sites across the transcriptome. These U2 complexes also contained the splicing regulators RBM5 and RBM10. We found RBM5 and RBM10 bound to nearly all branch site complexes and not simply those at regulated exons. The deletion of a conserved RBM5/RBM10 peptide sequence, including a zinc finger motif, disrupted U2 interaction and rendered the proteins inactive for the repression of many alternative exons. We propose a model where RBM5 and RBM10 regulate splicing as components of the U2 snRNP complex following branch site base pairing.

InPACT: a computational method for accurate characterization of intronic polyadenylation from RNA sequencing data.

Alternative polyadenylation can occur in introns, termed intronic polyadenylation (IPA), has been implicated in diverse biological processes and diseases, as it can produce noncoding transcripts or transcripts with truncated coding regions. However, a reliable method is required to accurately characterize IPA. Here, we propose a computational method called InPACT, which allows for the precise characterization of IPA from conventional RNA-seq data. InPACT successfully identifies numerous previously unannotated IPA transcripts in human cells, many of which are translated, as evidenced by ribosome profiling data. We have demonstrated that InPACT outperforms other methods in terms of IPA identification and quantification. Moreover, InPACT applied to monocyte activation reveals temporally coordinated IPA events. Further application on single-cell RNA-seq data of human fetal bone marrow reveals the expression of several IPA isoforms in a context-specific manner. Therefore, InPACT represents a powerful tool for the accurate characterization of IPA from RNA-seq data.

Transcription and splicing variations of SR genes accompany with genome-wide accumulation of long-introns in pine.

Most of mRNAs in Eukaryote were matured after the removal of introns in their pre-mRNA transcripts. Serine/arginine-rich (SR) proteins are a group of splicing regulators regulating the splicing processes globally. Expressions of SR proteins themselves were extensively regulated, at both transcription and splicing levels, under different environmental conditions, specially heat stress conditions. The pine genome is characterized by super-long and easily methylated introns in a large number of genes that derived from the extensive accumulation of transposons (TEs). Here, we identified and analyzed the phylogenetic characteristics of 24 SR proteins and their encoding genes from the pine genome. Then we explored transcription and pre-mRNA splicing expression patterns of SR genes in P. massoniana seedlings under normal and heat stress temperature conditions. Our results showed that the transcription patterns of SR genes in pine exhibited significant changes compared to other plant species, and these changes were not strictly correlated with the intron length and DNA methylation intensity of the SR genes. Interestingly, none of the long introns of SR genes underwent alternative splicing (AS) in our experiment. Furthermore, the intensity of AS regulation may be related to the potential DNA methylation intensity of SR genes. Taken together, this study explores for the first time the characteristics of significant variations in the transcription and splicing patterns of SR proteins in a plant species with an over-accumulation of super-long introns.

Distinct functions for the paralogous RBM41 and U11/U12-65K proteins in the minor spliceosome.

Here, we identify RBM41 as a novel unique protein component of the minor spliceosome. RBM41 has no previously recognized cellular function but has been identified as a paralog of U11/U12-65K, a known unique component of the U11/U12 di-snRNP. Both proteins use their highly similar C-terminal RRMs to bind to 3'-terminal stem-loops in U12 and U6atac snRNAs with comparable affinity. Our BioID data indicate that the unique N-terminal domain of RBM41 is necessary for its association with complexes containing DHX8, an RNA helicase, which in the major spliceosome drives the release of mature mRNA from the spliceosome. Consistently, we show that RBM41 associates with excised U12-type intron lariats, is present in the U12 mono-snRNP, and is enriched in Cajal bodies, together suggesting that RBM41 functions in the post-splicing steps of the minor spliceosome assembly/disassembly cycle. This contrasts with U11/U12-65K, which uses its N-terminal region to interact with U11 snRNP during intron recognition. Finally, while RBM41 knockout cells are viable, they show alterations in U12-type 3' splice site usage. Together, our results highlight the role of the 3'-terminal stem-loop of U12 snRNA as a dynamic binding platform for the U11/U12-65K and RBM41 proteins, which function at distinct stages of the assembly/disassembly cycle.

Understanding and Rescuing the Splicing Defect Caused by the Frequent ABCA4 Variant c.4253+43G>A Underlying Stargardt Disease.

Pathogenic variants in ABCA4 are the underlying molecular cause of Stargardt disease (STGD1), an autosomal recessive macular dystrophy characterized by a progressive loss of central vision. Among intronic ABCA4 variants, c.4253+43G>A is frequently detected in STGD1 cases and is classified as a hypomorphic allele, generally associated with late-onset cases. This variant was previously reported to alter splicing regulatory sequences, but the splicing outcome is not fully understood yet. In this study, we attempted to better understand its effect on splicing and to rescue the aberrant splicing via antisense oligonucleotides (AONs). Wild-type and c.4253+43G>A variant-harboring maxigene vectors revealed additional skipping events, which were not previously detected upon transfection in HEK293T cells. To restore exon inclusion, we designed a set of 27 AONs targeting either splicing silencer motifs or the variant region and screened these in maxigene-transfected HEK293T cells. Candidate AONs able to promote exon inclusion were selected for further testing in patient-derived photoreceptor precursor cells. Surprisingly, no robust splicing modulation was observed in this model system. Overall, this research helped to adequately characterize the splicing alteration caused by the c.4253+43G>A variant, although future development of AON-mediated exon inclusion therapy for ABCA4 is needed.

Cytochrome c oxidase subunit I gene in Thalassiosira nordenskioeldii strains inhabiting in cold and warm sea waters.

From the biota beneath the sea ice in Lake Saroma, which is adjacent to Sea of Okhotsk, a diatom culture of Saroma 16 was isolated. Strutted processes and a labiate process in Saroma 16 were characteristic of those in Thalassiosira nordenskioeldii. Similarity search analysis showed that the 826-bp rbcL-3P region sequence of this strain was 100% identical to multiple sequences registered as T. nordenskioeldii in a public database. The 4305-bp PCR-amplified mitochondrial cytochrome c oxidase subunit I (COI) gene (COI)-5P region of Saroma 16 included a 1060-bp open reading frame (ORF), which was interrupted by 934-bp and 2311-bp introns that included frame-shifted ORFs encoding reverse-transcriptase (RTase)-like proteins. Previous reports showed that a strain of the same species, CNS00052, originating from the East China Sea included no introns in the COI, whereas North Atlantic Ocean strains of the same species, such as CCMP992, CCMP993, and CCMP997, included a 2.3-kb intron in the same position as Saroma 16.

Research Progress of Group II Intron Splicing Factors in Land Plant Mitochondria.

Mitochondria are important organelles that provide energy for the life of cells. Group II introns are usually found in the mitochondrial genes of land plants. Correct splicing of group II introns is critical to mitochondrial gene expression, mitochondrial biological function, and plant growth and development. Ancestral group II introns are self-splicing ribozymes that can catalyze their own removal from pre-RNAs, while group II introns in land plant mitochondria went through degenerations in RNA structures, and thus they lost the ability to self-splice. Instead, splicing of these introns in the mitochondria of land plants is promoted by nuclear- and mitochondrial-encoded proteins. Many proteins involved in mitochondrial group II intron splicing have been characterized in land plants to date. Here, we present a summary of research progress on mitochondrial group II intron splicing in land plants, with a major focus on protein splicing factors and their probable functions on the splicing of mitochondrial group II introns.

Methylation Patterns of the FKBP5 Gene in Association with Childhood Maltreatment and Depressive Disorders.

Childhood maltreatment is an important risk factor for adult depression and has been associated with changes in the hypothalamic pituitary adrenal (HPA) axis, including cortisol secretion and methylation of the FKBP5 gene. Furthermore, associations between depression and HPA changes have been reported. This study investigated the associations of whole-blood FKBP5 mRNA levels, serum cortisol levels, childhood maltreatment, and depressive symptoms with the whole-blood methylation status (assessed via target bisulfite sequencing) of 105 CpGs at the FKBP5 locus using data from the general population-based Study of Health in Pomerania (SHIP) (N = 203). Both direct and interaction effects with the rs1360780 single-nucleotide polymorphism were investigated. Nominally significant associations of main effects on methylation of a single CpG site were observed at intron 3, intron 7, and the 3'-end of the gene. Additionally, methylation at two clusters at the 3'-end and intron 7 were nominally associated with childhood maltreatment × rs1360780 and depressive symptoms × rs1360780, respectively. The results add to the understanding of molecular mechanisms underlying the emergence of depression and could aid the development of personalised depression therapy and drug development.

Nutrient-dependent regulation of a stable intron modulates germline mitochondrial quality control.

Mitochondria are inherited exclusively from the mothers and are required for the proper development of embryos. Hence, germline mitochondrial quality is highly regulated during oogenesis to ensure oocyte viability. How nutrient availability influences germline mitochondrial quality control is unclear. Here we find that fasting leads to the accumulation of mitochondrial clumps and oogenesis arrest in Drosophila. Fasting induces the downregulation of the DIP1-Clueless pathway, leading to an increase in the expression of a stable intronic sequence RNA called sisR-1. Mechanistically, sisR-1 localizes to the mitochondrial clumps to inhibit the poly-ubiquitination of the outer mitochondrial protein Porin/VDAC1, thereby suppressing p62-mediated mitophagy. Alleviation of the fasting-induced high sisR-1 levels by either sisR-1 RNAi or refeeding leads to mitophagy, the resumption of oogenesis and an improvement in oocyte quality. Thus, our study provides a possible mechanism by which fasting can improve oocyte quality by modulating the mitochondrial quality control pathway. Of note, we uncover that the sisR-1 response also regulates mitochondrial clumping and oogenesis during protein deprivation, heat shock and aging, suggesting a broader role for this mechanism in germline mitochondrial quality control.

Chloroplast genome analyses of Caragana arborescens and Caragana opulens.

BACKGROUND: Numerous species within the genus Caragana have high ecological and medicinal value. However, species identification based on morphological characteristics is quite complicated in the genus. To address this issue, we analyzed complete plastid genome data for the genus. RESULTS: We obtained chloroplast genomes of two species, Caragana arborescens and Caragana opulens, using Illumina sequencing technology, with lengths of 129,473 bp and 132,815 bp, respectively. The absence of inverted repeat sequences in the two species indicated that they could be assigned to the inverted repeat-lacking clade (IRLC). The genomes included 111 distinct genes (4 rRNA genes, 31 tRNA genes, and 76 protein-coding genes). In addition, 16 genes containing introns were identified in the two genomes, the majority of which contained a single intron. Repeat analyses revealed 129 and 229 repeats in C. arborescens and C. opulens, respectively. C. arborescens and C. opulens genomes contained 277 and 265 simple sequence repeats, respectively. The two Caragana species exhibited similar codon usage patterns. rpl20-clpP, rps19-rpl2, and rpl23-ycf2 showed the highest nucleotide diversity (pi). In an analysis of sequence divergence, certain intergenic regions (matK-rbcL, psbM-petN, atpA-psbI, petA-psbL, psbE-petL, and rps7-rps12) were highly variable. A phylogenetic analysis showed that C. arborescens and C. opulens were related and clustered together with four other Caragana species. The genera Astragalus and Caragana were relatively closely related. CONCLUSIONS: The present study provides valuable information about the chloroplast genomes of C. arborescens and C. opulens and lays a foundation for future phylogenetic research and molecular marker development.

Genome Identification and Expression Profiling of the PIN-Formed Gene Family in Phoebe bournei under Abiotic Stresses.

PIN-formed (PIN) proteins-specific transcription factors that are widely distributed in plants-play a pivotal role in regulating polar auxin transport, thus influencing plant growth, development, and abiotic stress responses. Although the identification and functional validation of PIN genes have been extensively explored in various plant species, their understanding in woody plants-particularly the endangered species Phoebe bournei (Hemsl.) Yang-remains limited. P. bournei is an economically significant tree species that is endemic to southern China. For this study, we employed bioinformatics approaches to screen and identify 13 members of the PIN gene family in P. bournei. Through a phylogenetic analysis, we classified these genes into five sub-families: A, B, C, D, and E. Furthermore, we conducted a comprehensive analysis of the physicochemical properties, three-dimensional structures, conserved motifs, and gene structures of the PbPIN proteins. Our results demonstrate that all PbPIN genes consist of exons and introns, albeit with variations in their number and length, highlighting the conservation and evolutionary changes in PbPIN genes. The results of our collinearity analysis indicate that the expansion of the PbPIN gene family primarily occurred through segmental duplication. Additionally, by predicting cis-acting elements in their promoters, we inferred the potential involvement of PbPIN genes in plant hormone and abiotic stress responses. To investigate their expression patterns, we conducted a comprehensive expression profiling of PbPIN genes in different tissues. Notably, we observed differential expression levels of PbPINs across the various tissues. Moreover, we examined the expression profiles of five representative PbPIN genes under abiotic stress conditions, including heat, cold, salt, and drought stress. These experiments preliminarily verified their responsiveness and functional roles in mediating responses to abiotic stress. In summary, this study systematically analyzes the expression patterns of PIN genes and their response to abiotic stresses in P. bournei using whole-genome data. Our findings provide novel insights and valuable information for stress tolerance regulation in P. bournei. Moreover, the study offers significant contributions towards unraveling the functional characteristics of the PIN gene family.

On the genetic basis of tail-loss evolution in humans and apes.

The loss of the tail is among the most notable anatomical changes to have occurred along the evolutionary lineage leading to humans and to the 'anthropomorphous apes'1-3, with a proposed role in contributing to human bipedalism4-6. Yet, the genetic mechanism that facilitated tail-loss evolution in hominoids remains unknown. Here we present evidence that an individual insertion of an Alu element in the genome of the hominoid ancestor may have contributed to tail-loss evolution. We demonstrate that this Alu element-inserted into an intron of the TBXT gene7-9-pairs with a neighbouring ancestral Alu element encoded in the reverse genomic orientation and leads to a hominoid-specific alternative splicing event. To study the effect of this splicing event, we generated multiple mouse models that express both full-length and exon-skipped isoforms of Tbxt, mimicking the expression pattern of its hominoid orthologue TBXT. Mice expressing both Tbxt isoforms exhibit a complete absence of the tail or a shortened tail depending on the relative abundance of Tbxt isoforms expressed at the embryonic tail bud. These results support the notion that the exon-skipped transcript is sufficient to induce a tail-loss phenotype. Moreover, mice expressing the exon-skipped Tbxt isoform develop neural tube defects, a condition that affects approximately 1 in 1,000 neonates in humans10. Thus, tail-loss evolution may have been associated with an adaptive cost of the potential for neural tube defects, which continue to affect human health today.